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We develop the theory for shallow-core to valence excitations when the multiplet spread is larger 
than the core-hole attraction, e.g., if the core and valence orbitals have the same principal quantum 
number. This results in a cross-over from bound to virtual-bound excited states with increasing 
energy and in large differences between dipole and high-order multipole transitions, as observed in 
inelastic x-ray scattering. The theory is important to obtain ground state information from x-ray 
spectroscopies of strongly correlated transition metal, rare-earth and actinide systems. 
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The actinides and their compounds are attracting seri- 
ous attention from the condensed matter community due 
to their exotic properties []]]. Examples of these are the 
extremely rich phase diagram of Plutonium Q and the 
unsolved "hidden order" transition [3j in URu2Si2- Their 
properties interpolate between more itinerant 3d transi- 
tion metal (TM) systems and more localized Ad TM com- 
pounds 4], and exhibit strong interplay between spin, 
charge, orbital, and lattice degrees of freedom. Thus, ver- 
satile experimental techniques are needed to unravel the 
physics operative in these systems. Core-level spectro- 
scopies, like x-ray absorption spectroscopy (XAS), have 
been extremely successful in providing information re- 
garding the ground state of TM and rare-earth (RE) sys- 
tems, relying largely on theoretical interpretations based 
on local correlated models with full multiplet effects [B[ 
and atomic selection rules. 

The success of such local multiplet models relies 
strongly on two facts : (1) their electronic structure is 
largely governed by local correlation physics and point 
group symmetry, and (2) the final state core-hole strongly 
binds the extra d or f electron, so that all core-valence 
multiplets (CVM) form excitonic bound states. However, 
core-valence excitations within the same principal quan- 
tum number (n)-shell, like id-Af transitions in RE, or 
5d-5f transitions in the actinides, pose a problem be- 
cause the CVM spread is ~ 20-25 eV, often much larger 
than the average core-hole valence-electron attractive po- 
tential (Q) itself! This places higher lying terms up 
in the conduction band, 'autoionizing' the extra / elec- 
tron. This mixing with continua gives rise to very broad, 
virtual-bound (V-B) Fano-resonances Q, that cannot be 
interpreted in terms of local models. We note that this 
is an additional effect to those involving the decay of the 
core- hole itself, which have been elegantly described 0, IH 
for these giant dipole resonances (GDR) [l"oj |. 

The strong hybridization of the 5/ electrons with con- 
duction band states requires detailed information about 



interatomic interactions and band structure effects, rul- 
ing out the usefulness of atomistic models. In this pa- 
per, we emphasize that the structure of core-level exci- 
tations within the same n-shell, is fundamentally differ- 
ent from that between different n-shells, with the exam- 
ple of the 5d— >5f non-resonant inelastic x-ray scattering 
(NIXS) in the actinides. We also show that the high- 
multipole (HM) transitions in NIXS to strongly bound 
CVM states unlike the dipole restricted transitions 
in XAS and EELS, can still be treated within local mod- 
els, but with strongly renormalized parameters, due to 
very large configuration interaction (CI) in the final state. 

Unlike dipole restricted XAS, NIXS can access tran- 
sitions involving high-order multipoles as exemplified by 
the observation of d-d transitions in TM compounds [12|, 
13(. It also gives us more information of what the true 
ground state of the system actually was, even for dipole- 



allowed transitions [ll|, [14|. NIXS uses the first order 
scattering, off-resonance, due to the (e 2 /2mc 2 )A- A term 
in the light- matter coupling [HI . The corresponding dou- 
ble differential cross section is given by [l5| the Thompson 
scattering cross-section, times the material dependent dy- 
namical structure factor : 
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where, q—ki~kf is the photon momentum transfer, 
and hoj=Hu}i—Hujf is the energy loss. The transition 
operator can be multipole expanded [3] as : e l ^' r = 

where, ji(qr) are spherical Bessel functions, while 
Cm{9,(j)) are renormalized spherical harmonics, both of 
order I. Only terms with \lf — Zj| < I < (If + h) and 
(I + li + If) even, survive in the infinite sum. 

As a relevant illustration we show, in the inset-1 to 
Fig. 1(a), the radial transition probabilities for the vari- 
ous allowed channels in the bd^hf NIXS of the Th 4 + 
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FIG. 1: [color online) (a) The variation of the radial transi- 
tion probabilities with q, for the three component channels, 
I — l(dipole), 3(octupole) and 5(triakontadipole) (inset-1); 
plots of the 5d, 5/ and 6/ atomic H-F radial wavefunctions for 
Th 4+ (mset-2); and calculated NIXS (5(g,w)) for Th 4+ (5/°) 
including final-state CI with the 6/ level, for 60% of atomic 
C/(5/-6/) values (main). Both the 5/-like and 6/-like re- 
gions are shown. The bare component (1=1, 3, 5) spectra (not 
to scale) are shown at the bottom, (b) Schematic illustrating 
the high sensitivity of the dipole term vis-a-vis the HM terms, 
to the strongly term-dependent CI with the 6/. 

(5/°) system (relevant for TI1O2), plotted against q, 
as obtained using Cowan's atomic Hartree-Fock (H-F) 
code From the selection rules, this transition has 

allowed channels for the dipole (1=1), octupole (1=3) and 
the triakontadipole (1=5) sectors. The q and I depen- 
dence of the transition probabilities can be understood 
from the properties of ji(qr), as described by Haverkort 
et al. [H. 

While, details of the experimental data are discussed 



in the experimental paper 17J , here we briefly remind the 
reader of the key results obtained therein, viz. : (i) NIXS 
data for actinides show the dichotomy that HM features 
are sharp and excitonic, while the Fano-like asymmet- 
ric @ dipole feature is spread out over a very large en- 
ergy range; (ii) in stark contrast, the calculated low-q 



dipolar NIXS spectra are much sharper, intense and sim- 
ilar to the lower energy HM multiplets; (in) in order to 
obtain reasonable agreement with experimental spectra 
in terms of the peak positions, we needed to drastically 
scale down the 5d-5/ Coulomb and exchange Slater in- 
tegrals (F* f ,G k df ), to 60% (TI1O2) or 50% (U0 2 ) of their 
atomic H-F values! This is hard to justify for these rather 
atomic-like 5/ wavefunctions ^ n passing, we note 

that the dichotomy, discussed in (i) above, is also ob- 
served in the 4d^4/ NIXS of the RE and the 3p->3d 
NIXS of TM compounds [14J, and between the "doubly 
forbidden" pre-edge peak and the GDR, in the O45 XAS 
of actindes showing that this is a generic feature for 
shallow core-valence transitions within the same n-shell. 

To understand the origin of the apparent large reduc- 
tions in the Slater integrals, we plot in the inset-1 of Fig. 
1(a), the 5d, 5/ and the 6/ radial wavefunctions, ob- 
tained from an atomic H-F calculation for TI1O2 (Th 4+ , 
ground configuration 5/°). The 5d and 5/ orbitals (same 
n-shell) overlap very strongly, which accounts for the 
large values of the atomic (F^,G^) integrals. In con- 
trast, the more diffuse 6/ orbital, with an additional 
radial node, overlaps only weakly with the 5d. Hence, 
a final state (CI) between 5d 9 5f 1 and 5d 9 Qf 1 , via the 
C/(5/-6/) matrix elements [2(J, would effectively expand 
the radial part of the 5/ wavefunction, reducing the 5d- 
5/ Slater integrals. Now, the mixing depends also on the 
energy separation, AE, between the center-of-gravities 
(CG) of configurations involved, which is smaller in the 
actinides, than between 4d 9 4/ 1 and Ad 9 5f 1 in the RE. 
This results in a strong multiplet-dependent CI, because 
the multiplet splitting in hd?'z>f 1 is much larger than that 
in 5d 9 6f x , as shown schematically in Fig. 1(b). To illus- 
trate this point, the NIXS spectra for Th0 2 (Th 4+ ) are 
calculated (using the XTLS8.3 code [2lj) on the basis 
of the above model, as a coherent combination of the 
transitions hf-^bd^Sf 1 and 5/°^5d 9 6/ 1 [Hj]. Here we 
leave the (F^G^j) integrals unaltered at their atomic 
H-F values, while the scaling of CI(5f-6f) is varied to 
obtain agreement with experimental peak positions. This 
presents a more natural and physical mechanism for un- 
derstanding the reduced multiplet spread observed exper- 
imentally. The fact that we need to scale the CI integrals 
merely indicates that the band-like 6/ state is poorly 
approximated by the atomic H-F calculations, and that 
CI with numerous other states is neglected here. The 
full g-dependent spectra, at the optimized 60% scaling 
of C/(5/-6/)(Fig. 1(a)), is spread over a wide energy 
range and consists of 5 / -like and 6 / -like regions (marked 
in the figure), although most of the spectral weight lies 
in the 5/-like region due to dominance of radial matrix 
elements. The CI serves to reduce the 5d-5f multiplet 
spread considerably (while enhancing the 5d-6f spread), 
especially pushing the dipole peak close to the HM peaks, 
in good agreement with experiments 17]. Also the com- 
ponent spectra for the 1=1,3,5 channels plotted at the 
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FIG. 2: (color online) (a) Experimental g-averaged NIXS for 
TI1O2, compared with the component (1=1, 3, 5) summed cal- 
culated spectra for atomic values of (F%,G%), and the CI(5f- 
6/) varied from a scaling of 20% to their full atomic val- 
ues. Only peak-positions are relevant. The best agreement 
is obtained at 60% scaling of C7(5/-6/). Variation of the 
high-energy dipole peak position with varying scaling of : (b) 
C7(5/-6/), and (c) (F% f ,G k df ). While the latter shows a sim- 
ple linear trend, the former shows a parabolic behavior with 
signs of saturation. 

base of Fig. 1(a) (vertical scale not to be compared with 
the actual spectra above) clearly show that it also repro- 
duces the mild shift of the HM peak at ~ 93 eV, due to 
a g-dependent weight transfer between the 1=3 and 1=5 
channels, as seen in experiments 17]. The dipole states, 
in general, mix much more and respond much more sen- 
sitively to the CI than the HM states. This also means 
that there is a substantial amount of interference between 
the 5/-like and 6/-like dipole transitions. 

The origin of the higher sensitivity of the dipole com- 
pared to the HM is explained schematically in Fig. 1(b). 
Before CI (left), the difference (AE) in the CG ener- 
gies for the two multiplets, 5d 9 5f 1 (width W mu it) and 
hd^Qf 1 (width Wmuit), is ~ 23-24 eV, from H-F calcu- 
lations. Although both the multiplets involve exactly 
the same terms, W mu it (~ 25 eV) is much larger than 
Wmuit (~ 7 eV), demonstrating the difference between 
states involving the same versus different principle quan- 
tum numbers. Also, for a less-than-half-filled system with 
a low J ground state ( J=0 for TI1O2) , the highest multi- 
plets (red) are generally dipole- allowed, while the lowest 
ones (blue) are the HM-allowed terms. Now C/(5/-6/) 
only mixes terms of the same symmetry (l6j . e.g., the 



red (blue) states at the top (bottom) of hd^hf 1 , mix only 
with the red (blue) states at the top (bottom) of 5d 9 6/ x . 
The result after mixing is shown in the right panel. Since 
Wmuit « W mu it, the effective energy denominator for 
mixing of the red (dipole) terms is much smaller than 
that for the blue (HM) terms, causing the observed dif- 
ferences in the shifts. Due to these strong correlation 
effects, the effective "screening" becomes highly term de- 
pendent, explaining the strong asymmetry in the behav- 
ior of the dipole and the HM terms, which is not captured 
by a uniform reduction of the Slater integrals fl7i j . 

To show that the two approaches are qualitatively dif- 
ferent, we compare in Fig. 2(a) the experimental q- 
averaged NIXS (topmost), with the calculated sum of 
the three component spectra (1=1,3 and 5), keeping 
(F^,G%) fixed at their atomic values, while the CJ(5/- 
6/) are switched on and gradually increased to their 
atomic value (100%) (top to bottom). For this purpose 
only the peak positions are relevant. As already noted, 
a good agreement with HM peak positions is obtained 
for 60% reduction of CI(5f-6f). But more importantly, 
with gradual uniform reduction in (F%,G%) we would 
expect a linear movement of the dipole towards the HM 
features. On the other hand with changing degree of 
CI, the dipole moves towards the HM peaks in a non- 
linear manner, showing signs of saturation, as governed 
by level-repulsion physics. This contrasting behavior is 
shown in Figs. 2(b) and 2(c), where we have plotted the 
dipole peak position in the two cases, as a function of the 
scaling of (F^,G^) and of C/(5/-6/), respectively. 

We now turn to the problem of the experimentally 
observed low relative amplitude, large width and non- 
lorentzian line shape of the GDR. This can be under- 
stood on the basis of Figs. 3(a)-(b). Fig. 3(a) shows 
the NIXS final state, 5d 9 5/ n+1 , of a 5/™ actinide sys- 
tem. In the absence of the 5d core-hole, this would be 
identical to the inverse photoemission final state 5/™ +1 
(dashed line), and lies just below or within a continuum 
(shaded) . Primarily, the core-hole provides an attractive 
scalar potential, Q, which is often large enough to pull 
down the 5d 9 5/ n+1 state, out of the continuum, forming 
a bound core-hole exciton. However, in a more complete 
picture (Fig. 3(b)), the core-hole /-electron multipolc 
interaction, also yields a very broad multiplet structure 
(width Wmuit) about the CG of 5d 9 5/' i+1 . The Slater 
integrals are very large for these transitions occurring 
within the same n-shell, implying W m uit»Q- Thus 
while the HM states towards the bottom of the multiplet, 
still form bound states, the high-lying dipolar terms are 
pushed up into the continuum, offsetting the effect of Q, 
and their mixing gives rise to the GDR with characteris- 
tic Fano lineshapes, as seen experimentally. This physi- 
cal picture clearly demonstrates why the dipole and the 
HM states show a crossover from V-B to bound charac- 
ter, in the 5d->5/ (actinides) [HI, the 4(i— s-4/ (RE com- 
pounds) [ll[ , or the 3p— s-3d (TM compounds) edges [3| • 
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FIG. 3: [color online) Schematics showing (a) the effect of 
a scaler core-hole potential (Q) on the NIXS final state; and 
(b)its fanning out due to strong core- valence multiplet effects 
leading to a gradual change from bound-states (lower terms) 
to V-B resonances (higher terms) as a function of energy; (c) 
the S(q, cu) from a model calculation for Th 4+ (lines) including 
a Fano-effect of the mixing of the atomic 5d-5/ transition with 
that to a fictitious broad band, causing a broad, asymmetric 
lineshape of the high energy dipole term, that agrees well with 
experimental (/-dependent NIXS for TI1O2 Q3] (symbols). 



In Fig. 3(c) we show S(q,u) from a model calcula- 
tion for Th 4+ (lines) that includes transitions from the 
5d core-level to both the (5/, 6/) levels as before (al- 
lowed 1=1,3,5), and to a fictitious 7p-like discretized 
band (allowed Z=l, 3) . The band is so positioned that 
it starts below the high energy dipole state, but above 
the HM states. The 2=1,3 channels can interfere via 
the CI(5f-7p) matrix elements [H We find that 

while the lower lying 1= 3, 5 peaks remain sharp and 



excitonic, the high-energy dipole feature forms a GDR, 
just as discussed above. Interestingly, a dipole-allowed 
peak present at lower energy (~ 88 eV) is not broadened 
by this mechanism, implying that the position within the 
multiplet, rather than the symmetry of the state, decides 
its fate. A fairly good comparison with the experimental 
q-dependent NIXS for Th0 2 (symbols) is obtained if 
we use a somewhat larger Lorentzian width for the GDR 



than the HM states, in order to simulate the multiplet 
dependent core-hole decay probabilities, not included in 
the present calculation. It is important to note that this 
dichotomy, between the dipole and the HM, could be re- 
versed in cases where the dipole allowed states are lower 
in energy than the HM, like for more-than-half-fillcd sys- 
tems which have large ground state J values 22j | . 

In conclusion, the modeling of same n-shell NIXS is 
complicated by the simultaneous presence of V-B and 
bound states within the same final-state multiplet. The 
complex V-B resonances, involving non-local effects, pro- 
vide insight into the hybridization of the locally ex- 
cited core-electron with continua, and about core-hole 
decay processes. The dipole-forbidden bound states (not 
prominent in XAS) are modeled using a local, atomic 
CI approach, that provides direct ground state informa- 
tion. It also explains the apparent strong reduction of the 
atomic Slater integrals, which effectively become term 
dependent. The HM transitions and especially their an- 
gular dependence in single crystal studies [TJ, LL3| are 
also expected to provide detailed information on the im- 
portance of "orbital-ordering" in the ground state and 
changes at phase transitions in "hidden order" materi- 
als like URu2Si2 0]. The V-B resonances are modeled 
with an additional term, mixing the 5/ states with a 
conduction band continuum, resulting in their broad and 
Fano-like line shapes. More realistic approaches to the 
latter would include an energy-dependent hybridization 
with a realistic density of states, and explicit core-hole 
decay processes. This is a topic of future investigations. 
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